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We study a random even subgraph of a finite graph G with a general edge- weight 
p £ (0, 1). We demonstrate how it may be obtained from a certain random-cluster 
measure on G, and we propose a sampling algorithm based on coupling from the 
I past. A random even subgraph of a planar lattice undergoes a phase transition at 

' the parameter- value ^pc, where pc is the critical point of the q = 2 random-cluster 

model on the dual lattice. The properties of such a graph are discussed, and are 
^ i related to Schramm-Lowner evolutions (SLE). 



1 Introduction 

Our purpose in this paper is to study a random even subgraph of a finite graph G = (V, E), 
and to show how to sample such a subgraph. A subset F of is called even if, for all 
X & V, X is incident to an even number of elements of F. We call the subgraph (V, F) 
even if F is even, and we write S for the set of all even subsets F of E. It is standard that 
every even set F may be decomposed as an edge-disjoint union of cycles. Let p G [0, 1). 
The random even subgraph of G with parameter p is that with law 

Pp(^) = ^pl^'(l-P)l^\^', Fes, (Ll) 



THE ELECTRONIC JOURNAL OF COMBINATORICS 16 (2009), #R00 



1 



where — Zq{p) is the appropriate normahzing constant. 

We may express pp as follows in terms of product measure on E. Let (j)p be product 
measure with density p on the configuration space Q = {0, 1}^. For uj E and e E E, 
we call e uj-open if a;(e) = 1, and u-closed otherwise. Let duj denote the set of vertices 
X &V that are incident to an odd number of a;-open edges. Then 

where up is the edge-configuration whose open edge-set is F. In other words, (f)p describes 
the random subgraph of G obtained by randomly and independently deleting each edge 
with probability 1 — p, and Pp is the law of this random subgraph conditioned on being 
even. 

Random even graphs are closely related to the Ising model and the random-cluster 
model on G, and we review these models briefly. Let /? e (0, oo) and 

_2/3 2tanh/3 

p = 1 — e ^ — T-—. (1.3) 

^ l + tanh/3 ^ ' 

The Ising model on G has configuration space E = {—1, +1}^, and probability measure 

7r;3((7) = -^exp|/3^(7^(Ty|, creE, (1.4) 

where — Zq{(5) is the partition function that makes tt^j a probability measure, and 

e = {x.y) denotes an edge with endpoints y. A spin-duster of a configuration o" G E 
is a maximal connected subgraph of G each of whose vertices v has the same spin-value 
ay. A spin-cluster is termed a k cluster ii = k for all v belonging to the cluster. An 
important quantity associated with the Ising model is the 'two-point correlation function' 

y) = '^fi{(^x = (^y)- ^ = \T^I}{.(^x<^y), X,y eV, (1.5) 

where P{f) denotes the expectation of a random variable / under the probability measure 
P. 

The random-cluster measure on G with parameters p G (0, 1) and g = 2 is given as 
follows [it may be defined for general q > but we are concerned here only with the case 
q = 2]. Let 



= -^pl''('^)l(l -p)l^\''('")l2*^('^), cjen, (1.6) 
z 

where k{uj) denotes the number of u;-open components on the vertex-set V, i]{uj) = {e G 
£^ : a;(e) = 1} is the set of open edges, and Z^^ = Zq^{p) is the appropriate normalizing 
factor. 
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The relationship between the Ising and random-cluster models on G is well established, 
and hinges on the fact that, in the notation introduced above. 



y) 



where {x ^ y} is the event that x and y are connected by an open path. See [12] for 
an account of the random-cluster model. There is a relationship between the Ising model 
and the random even graph also, known misleadingly as the 'high-temperature expansion'. 
This may be stated as follows. For completeness, we include a proof of this standard fact 
at the end of the section, see also [3]. 

Theorem 1.7. Let 2p = 1 — e'"^^ where p G (0, |), and consider the Ising model with 
inverse temperature (3. Then 



A corresponding conclusion is valid for the product of cTx^ over any even family of 
distinct Xi e V . 

This note is laid out in the following way. In Section [2] we define a random even 
subgraph of a finite or infinite graph, and we explain how to sample a uniform even 
subgraph. In Section [3] we explain how to sample a non-uniform random even graph, 
starting with a sample from a random-cluster measure. An algorithm for exact sampling 
is presented in Section H] based on the method of coupling from the past. The structure 
of random even subgraphs of the square and hexagonal lattices is summarized in Section 



In a second paper [H], we study the asymptotic properties of a random even subgraph 
of the complete graph Kn- Whereas the special relationship with the random-cluster and 
Ising models is the main feature of the current work, the analysis of [H] is more analytic, 
and extends to random graphs whose vertex degrees are constrained to lie in any given 
subsequence of the non-negative integers. 

Remark 1.8. The definition (11. ip may be generalized by replacing the single parameter 
p by a family p = (pe : e G -E), just as sometimes is done for the random-cluster measure 
( II. 6p . see for example [20]; we let 





El 





For simplicity we will mostly consider the case of a single p. 



Proof of Theorem \1.7\ For cr G S, G ^2, let 
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where deg{v,uj) is the degree of v in the 'open' graph (y,ri{uj)). Then 



0)60 e=(i;,ui) e=(v,w) 

= e-^l^l exp I /3 J] a^a^ , a G S. (1.11) 

y e={v,w) j 

Similarly, 

Y,Z,{a,uj) = 2l^lpl''(-)l(l -p)'^\^(-)'l{a.=^}, ouen, (1.12) 

and 



J]fr,a,Zp(a,^) = 2l^lpl''(-)l(l - p)l^\''(-)ll{a^=|,,,|}, cu e n. (1.13) 



o-es 



By dniD, 

7r/3,2(0-;rO-s/j ^ 



and the claim follows by ffrT2D - ffrT3D . □ 



2 Uniform random even subgraphs 
2.1 Finite graphs 

In the case p = | in (11. ip . every even subgraph has the same probability, so pi describes 
a uniform random even subgraph of G. Such a random subgraph can be oljtained as 
follows. 

We identify the family of all spanning subgraphs of G = (V, E) with the family 2^ 
of all subsets of E. This family can further be identified with {0,1}^ = Z|^, and is 
thus a vector space over Z2; the addition is componentwise addition modulo 2 in {0, 1}^, 
which translates into taking the symmetric difference of edge-sets: Fi + F2 = Fi A F2 for 
Fi,F2 C E. 

The family of even subgraphs of G forms a subspace S of this vector space {0, 1}^, 
since Fi + F2 = Fi A F2 is even if Fi and F2 are even. (In fact, S is the cycle space 
Zi in the Z2-homology of G as a simplicial complex.) In particular, the number of even 
subgraphs of G equals 2^^^"^^ where c{G) = dim(£^); c{G) is thus the number of independent 
cycles in G, and, as is well known, 

c{G) = \E\-\V\ + k{G). (2.1) 
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Proposition 2.2. Let Ci,...,Cc be a maximal set of independent cycles in G. Let 
^i,...,^c be independent Be(|) random variables (i.e., the results of fair coin tosses). 
Then is a uniform random even subgraph of G. 

Proof. Ci, Cc is a basis of tlie vector space £ over Z2. □ 

One standard way of clioosing Ci, . . . , Cc is exploited in the next proposition. Another, 
for planar graphs, is given by the boundaries of the finite faces; this will be used in Section 
[5l In the following proposition, we use the term spanning subforest of G to mean a maximal 
forest of G, that is, the union of a spanning tree from each component of G. 

Proposition 2.3. Let (V, F) be a spanning subforest of G. Each subset X of E\F can be 
completed by a unique Y C F to an even edge-set Ex = X UY E S . Choosing a uniform 
random subset X ^ E \ F thus gives a uniform random even subgraph Ex of G. 

Proof. It is easy to see, and well known, that each edge Ci E E\F can be completed by 
edges in F to a unique cycle Gi\ these cycles form a basis of £ and the result follows by 
Proposition 12.21 (It is also easy to give a direct proof.) □ 

2.2 Infinite graphs 

Here, and only here, we consider even subgraphs of infinite graphs. Let G = {V, E) be 
a locally finite, infinite graph. We call a set c 2^ finitary if each edge in E belongs 
to only a finite number of elements in JF. If G is countable (for example, if G is locally 
finite and connected), then any finitary JF is necessarily countable. If JF c 2^ is finitary, 
then the (generally infinite) sum Yl^^-pX is a well-defined element of 2^, by considering 
one coordinate (edge) at a time; if, for simplicity, JF = {x, : i G /}, then includes 
a given edge e if and only if e lies in an odd number of the Xi. 

We can define the even subspace S of 2^ as before. (Note that we need G to be locally 
finite in order to do so.) If JF is a finitary subset of S, then Ylxey^-'^ ^ ^■ 

A finitary basis of £^ is a finitary subset J-" G S such that every element of S is the sum 
of a unique subset JF' C JF; in other words, if the hnear (over Z2) map 2-^ ^ £ defined by 
summation is an isomorphism. (A finitary basis is not a vector-space basis in the usual 
algebraic sense since the summations are generally infinite.) 

We define an infinite cycle in G to be a subgraph isomorphic to Z, i.e., a doubly infinite 
path. (It is natural to regard such a path as a cycle passing through infinity.) Note that, 
if F is an even subgraph of G, then every edge e E F belongs to some finite or infinite 
cycle in F: if no finite cycle contains e, removal of e would disconnect the component of 
F that contains e into two parts; since F is even both parts have to be infinite, so there 
exist infinite rays from the endpoints of e, which together with e form an infinite cycle. 

Proposition 2.4. The space £ has a finitary basis. We may choose such a finitary basis 
containing only finite or infinite cycles. 
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Proof. It suffices to consider the case when G is connected, and hence countable. We 
construct a finitary basis by induction. Order the edges in a fixed but arbitrary way as 
ei < 62 < ■ ■ ■ . Let hi be the first edge that belongs to an even subgraph of G, and 
choose a (finite or infinite) cycle Ci containing hi. Having chosen hi, Ci, . . . , Cn, 
consider the subspace Sn of all even subgraphs of G containing none of hi, ... , hn. If 
Sn = {0}, we stop, and write JF = {Ci, C2, . . . , Cn}- Otherwise, let hn+i be the earliest 
edge belonging to some non-trivial even subgraph Fn E £n, and choose a cycle Cn+i C F„ 
containing Either this process stops after finitely many steps, with the cycle set 

J-', or it continues forever, and we write for the countable set of cycles thus obtained. 
Finally, write H = {hi, h2, . . .}. We shall assume that H ^ 0, since the proposition is 
trivial otherwise. 

We claim that JF is a finitary basis for S. Note that 

hn G Gn, hj ^ Gn foT j < U. (2.5) 

Let e G -E, say e = e.^. If = hg for some s, then lies in only finitely many of the 
Gj. \{ Cf E E \ H and hg < < /i^+i for some s (or hg < e^. for all s), then lies in no 
member of Es-, so that it lies in only finitely many of the Gj. If Cr < hi, then lies in no 
Gj. In conclusion, is finitary. 

Next we show that no element F E S has more than one representation in terms 
of JF. Suppose, on the contrary, that ^j^jCj = '^^ipiGi. Then the sum of these two 
summations is the empty set. By (12.51) . there is no non-trivial linear combination of the 
Gi that equals the empty set, and therefore C,i = ipi for every i. 

Finally, we show that JF spans S, which is to say that the map 2-^ S defined by 
summation has range S. Let JF be the subspace of S spanned by JF. For H' C H, there 
is a unique element F' E T such that F' Ci H = H'; F' is obtained by an inductive 
construction that considers the Gj in order of increasing j, and includes a given Gj if: 
either hj G H' and hj lies in an even number of the Gi already included, or hj ^ H' and 
hj lies in an odd number of the Gi already included. 

Let F E S. By the above, there is a unique element F' E T satisfying F'dH = FCiH. 
Thus, F + F' is an even subgraph having empty intersection with H. Let be the earliest 
edge in F + F', if such an edge exists. Since G F + F', there exists s with hg < e^. 
With s chosen to be maximal with this property, we have that lies in no even subgraph 
of Sg, in contradiction of the properties of F + F' . Therefore, no such Cr exists, so that 
F + F' = 0, and F = F' eT as required. □ 

Given any finitary basis JF = {Ci, G2, . . .} of S, we may sample a uniform random 
even subgraph of G by extending the recipe of Proposition 2.2 to infinite sums: we let 
^1,^2, • • • be independent Be(|) random variables and take X^i'CjC'j- other words, we 
take the sum of a random subset of the finitary basis JF obtained by selecting elements 
independently with probability ^ each. Denote by p the ensuing probability measure on 
£. 

It turns out that p is specified in a natural way by its projections. Let Ei be a finite 
subset of E. The natural projection tce^ : {0, 1}^ — > {0, l}'^^ given by iiEiiuj) = (co'e)eG£;i 
maps £ onto a subspace £ei = t^Ei{£) of {0, 1}^^ 
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Theorem 2.6. Let G be a locally finite, infinite graph. The measure p given above is the 
unique probability measure on Q = {0, 1}^ such that, for every finite set Ei C E with 
£e^ 7^ 0, (co'e)egEi is Uniformly distributed on i-^-, 

pin^liA)) = IAhSeA/ISeA, AC {0,1}^^ (2.7) 

Proof. We may assume that G is connected since, if not, any p satisfying (12. 7p is a 
product measure over the different components of G. Note that every connected, locaUy 
finite graph is countable. 

We show next that there is a unique probabihty measure satisfying (12. 7p . This equation 
specifies its value on any cylinder event. By the Kolmogorov extension theorem, it suffices 
to show that this specification is consistent as Ei varies, which amounts to showing that 
a El C E2 C E with Ei,E2 finite, then the projection ttejEi '■ {0, 1}^^ {0, 1}^^ maps 
the uniform distribution on Se2 to the uniform distribution on Se^ ■ This is an immediate 
consequence of the fact that nE2Ei is a linear map of £e2 onto Se^- 

Finally we show that p satisfies (12.71) . Let Ei G E he finite. Since is finitary, its 
subset J-'i, containing cycles that intersect Ei, is finite. Since p is obtained from uniform 
product measure on JF, its projection onto Ei is uniform (on its range) also. □ 

Diestel [3, Chap. 8] discusses related results for the space of subgraphs spanned by the 
finite cycles, and relates them to closed curves in the Freudenthal compactification of G 
obtained by adding ends to the graph. It is tempting to guess that there may be similar 
results for even subgraphs and the one-point compactification of G (where all ends are 
identified to a single point at infinity). We do not explore this here, except to note that 
the finite and infinite cycles are exactly those subsets of the one-point compactification 
that are homeomorphic to a circle. 

3 Random even subgraphs via coupling 

We return to the random even subgraph with parameter p G [0, 1) defined by (II. ip for a 
finite graph G = {V,E). We show next how to couple the q = 2 random-cluster model 
and the random even subgraph of G. Let p G [0, |], and let be a realization of the 
random-cluster model on G with parameters 2p and q = 2. Let R = (V, 7) be a uniform 
random even subgraph of {V,t]{uj)). 

Theorem 3.1. Let p G [0, |]. The graph R = (V, 7) is a random even subgraph of G with 
parameter p. 

This recipe for random even subgraphs provides a neat method for their simulation, 
provided p < ^. One may sample from the random-cluster measure by the method of 
coupling from the past (see [21] and Section H]), and then sample a uniform random even 
subgraph by either Proposition 12.21 or Proposition 12.31 
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Proof. Let g C E he even. By the observations in Section [2.11 with c{uj) = c(V,?7(ct;)) 
denoting the number of independent cycles in the open subgraph, 

if g C r]{u), 

otherwise, 



= 9 \ ^) 
so that 



n7 = 9)= Yl 2-^^"V2p,2M. 

u}:gCr]{uj) 

Now c{uj) = \r]{uj) \ -\V\ + k{uj), so that, by (TO . 

oc ^ pl'?HI(i_2p)I^V';HI 

<^:gCj7(a;) 

= [p+ (1 -2p)]l^\^'lpl^'l 
= pl3l(l-p)l^\»l, ^CE. 
The claim follows. □ 

Let p E (|, !)• If G is even, we can sample from pp by first sampling a subgraph {V, F) 
from pi_p and then taking the complement {V,E \ F), which has the distribution pp. If 
G is not even, we adapt this recipe as follows. For W V and H O E, we say that H 
is W-even if each component of {V,H) contains an even number of members of W. Let 
W ^ he the set of vertices of G with odd degree, so that, in particular, E is ly-even. 
Let = {uj E Q : ri{uj) is W-even}. For u G , we pick disjoint subsets P* = P^, 
2 = 1, 2, ... , ^\W\, of r7(co'), each of which constitutes an open non-self-intersecting path 
with distinct endpoints lying in W, and such that every member of W is the endpoint of 
exactly one such path. Write P^j = [JiP^- 

Let r = 2(1 — p), and let (f)^2 be the random-cluster measure on Q with parameters 
r and q = 2 conditional on the event fl^ . We sample from 0^ to obtain a subgraph 
{V, ri{u))), from which we select a uniform random even subgraph {V, 7) by the procedure 
of the previous section. 

Theorem 3.2. Letp G (|, 1). T/ie graph S = {V, -E'\(7 A P^)) is a random even subgraph 
of G with parameter p. 

The recipes in Theorems 13.11 and 13.21 can be combined as follows. Consider the gen- 
eralized model mentioned in Remark 11.81 with one parameter p^, G (0, 1) for each edge 
e E E. Let A = {e E E : Pf. > |}. Define r^. = 2pf. when e ^ A and r^. = 2(1 — p^) 
when e E A. (Thus < rg < 1.) Let W = Wa he the set of vertices that are A-odd, i.e., 
endpoints of an odd number of edges in A. Sample u from the random-cluster measure 
with parameters r = {r^ : e E E) and q = 2, conditioned on rjlu) being H^-even, let P^ he 
as above (for W = Wa), and sample a uniform random even subgraph (V, 7) of (V, r]{uj)). 
For a discussion of relevant sampling techniques, see Section HI 
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Theorem 3.3. The graph S = (y,'j A P^^ A A) is a random even subgraph of G with 
the distribution pp given in (11.91) . 

Note that Theorems 13.11 and 13.21 are special cases of Theorem 13.31 with A = and 
A = E respectively. We find it more illuminating to present the proof of Theorem 13.21 in 
this more general setup. 

Proof of Theorem \3.3l and thus of Theorem \3.2[ Let F = A P^ A A he the resulting 
edge-set, and note that 

r]{iu) D-f A P^ = F A A. (3.4) 

Furthermore, if F is even, then F A A has odd degree exactly at vertices in W = Wa'-, 
hence fl3.4l) implies that necessarily uo G 

Given an even edge-set f ^ E, we thus obtain F = / if we first choose uj G with 
^(^) ^ / A A and then (having chosen P^) select 7 as the even subgraph f A A A P^^. 
Hence, for every u G fl^ with r]{uj) ^ f A A, we have F{F = f \ u) = 2^^^"^^ and 
summing over such uj we find 

P(F = /) CX Yl 2-^("Vr,2(^) 

OC 2-"(^)2'^(") Yl rf'\l - re)^-'^(") 

OC 2-l''(")l JJr^(")(l -re) 

- E n(|)"'"(i-^=)-> 

- n (?) n (-?). 

eG/AA eifAA 

With le denoting the indicator function of the event {e G /}, this can be rewritten as 
P(F = /) oc \{{rj2Y^{l - rj2Y-'^ \{{rj2f-^^{l - rj2f^ 

e^A e&A 

e^A eeA 

eGB 
OC Pp(/). 

The claim follows. □ 

There is a converse to Theorem 13. 1[ Take a random even subgraph {V, F) oiG = {V, E) 
with parameter p < |. To each e ^ F, we assign an independent random colour, blue 
with probability p/ (1 —p) and red otherwise. Let H be obtained from F by adding in all 
blue edges. 
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Theorem 3.5. The graph (V, H) has law 
Proof. For h (1 E, 



N |J| / \ \h\J\ /. r, \ 

p \ I p \ I 1 — 2p^ 



JCh, J even 

ocpl'^l(l-2p)l^\'^lAr(/i), 



1 — pj \ 1 — p 



where N{h) is the number of even subgraphs of {V,h). As in the above proof, N{h) = 
2W-\^\+''W where k{h) is the number of components of {V, h), and the proof is complete. 

□ 

An edge e of a graph is called cyclic if it belongs to some cycle of the graph. 

Corollary 3.6. For p E [0, |] and e E E, 

Pp{e is open) = ^02p,2(e is a cyclic edge of the open graph). 

By summing over e E E, we deduce that the mean number of open edges under pp is 
one half of the mean number of cyclic edges under 02p,2- 

Proof. Let uj E Q and let C be a maximal family of independent cycles of a;. Let R = {V, 7) 
be a uniform random even subgraph of {V,r]{uj)), constructed using Proposition 12.21 and 
C. For e E E, let Mg be the number of elements of C that include e. If Me > 1, the 
number of these cycles of 7 that are selected in the construction of 7 is equally likely 
to be even as odd. Therefore, 



P(e e 7 I u) 



I ifMe>l, 



if Me = 0. 

The claim follows by Theorem 13.11 □ 



4 Sampling an even subgraph 

It was remarked earlier that Theorem 13. II gives a neat way of sampling an even subgraph 
of G according to the probability measure rjp with p < |. Simply use coupling-from-the- 
past (cftp) to sample from the random-cluster measure 02p,2, and then flip a fair coin once 
for each member of some maximal independent set of cycles of G. 

The theory of cftp was enunciated in ^21J and has received much attention since. We 
recall that an implementation of cftp runs for a random length of time T whose tail 
is bounded above by a geometric distribution; it terminates with probability 1 with an 
exact sample from the target distribution. The random-cluster measure is one of the 
main examples treated in |2T] . We do not address questions of complexity and runtime in 
the current paper, but we remind the reader of the discussion in ^21j of the relationship 
between the mean runtime of cftp to that of the underlying Gibbs sampler. 
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The situation is shghtly more comphcated when p > | and G is not itself even, since 
the conditioned random-cluster measure used in Theorems 13.21 and l3.3l is neither monotone 
nor anti-monotone. We indicate briefly in this section how to adapt the technique of cftp 
to such a situation. 

Let ii^ be a non-empty finite set, and let /i be a probability measure on the product 
space Vt = {0, l}"^. We call /i monotone (respectively, anti-monotone) if /i(le | ^e) is 
non-decreasing (respectively, non-increasing) in ^ e ^2. Here, le is the indicator function 
that e is open, and C,e is the configuration obtained from ^ on ii^ \ {e}. For e & E, ijj & Q, 
and 6 = 0, 1, we write ipe configuration that agrees with ip oS e and takes the value 

b on e. 

It is standard that cftp may be used to sample from /i if /i is monotone (see [HI [28]). 
and it is explained in p!3] how to adapt this when /x is anti-monotone. We propose 
below a mechanism that results in an exact sample from n without any assumption of 
(anti-)monotonicity. This mechanism may be cast in the more general framework of the 
'bounding chain' of |18], but, unlike in that work, it makes use of the fact that fl is 
partially ordered. A similar approach was proposed in [19] under the title 'dominated 
CFTP', in the context of the simulation of point processes. 

Write = {u E fl : ^{u) > 0}, the subset of fl on which /i is strictly positive, and 
assume for simplicity that S*^ is increasing, and that 1 G 5^, where 1 (respectively, 0) 
denotes the configuration of 'all 1' (respectively, 'all 0'). This assumption is valid in the 
current setting, but is not necessary for all that follows. 

We start with the usual Gibbs sampler for fi. This is a discrete-time Markov chain 
G = {Gn : n > 0) on the state space that updates as follows. Suppose Gn = i- A 
uniformly distributed member of E is chosen, e say, and also a random variable U with 
the uniform distribution on [0, 1]. Then Gn+i = C where df) = ^(/) for f ^ e, and 

^'(e) = l^ if?7>/i(le|ee), 
\l iiU<fi{U\Q. 

The transition rule is well defined whenever G S*^. It is convenient to use the device of 
[T5] to extend this definition to configurations not in S^, and to this end we set 

/i(le I Q = max{/i(le I Ipe) ■ > ^e, i'l ^ S^} (4.1) 

when ^ S^. There is a degree of arbitrariness about this definition, which we follow 
for consistency with [T5] . 

Let (en, Un) be an independent sequence as above. Let {An, -B„ : n > 0) be a Markov 
chain with state space fi^, and {Ao,Bo) = (0, 1). Suppose (y4„,i?„) = {C,,ri) where ^ < //. 
We set {An+u i?„+i) = (C, v') where Cif) = e(/), v'if) = vif) for / ^ e„+i. At e = e„+i 
we set 

,^'(e) = 1 if and only if f/„+i < a, 
rj'{e) = 1 if and only if t/„+i < /3, 
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where 

a = a{^,T]) = min{/i(le | tpe) ■ < < Ve}, 2) 

P = Pi^^V) = max{/i(le I Ipe) ■■^e<i>e< Ve} ■ 

Since a < P, we have that ^' < t]'. 

We run the chain {A, B) starting at negative times, in the manner prescribed by cftp, 
and let T be the coalescence time. More precisely, for m > 0, let {Ak{m),B}~{m) : —m < 
< 0) denote the chain beginning with A_m{jn) = 0, B_^{m) = 1, using a fixed random 
sequence (e„, ?7n)° oo set 

T = min{m > : v4o(m) = Bo{m)}, 

so that v4o(T) = Bo{T). 

Theorem 4.3. If is increasing and 1 G S*^, then P{T < 00) = 1, and Aq{T) has law 

Proof. We prove only that P{T < 00) = 1. The second part is a standard exercise in 
cftp, and is easily derived as in [151 Thm 2.2]. By the definition of S"^ and (14. ip . there 
exists rj = ri{E,fi) > such that yu(le \ ^e) ^ V for all e G -E and ^ G In any given 
time-interval of length \E\, there is a strictly positive probability that the corresponding 
sequence (cj, Ui) satisfies E = {ci} and Ui < rj for all i. On this event, the lower process 
A takes the value 1 after the interval is past, so that coalescence has taken place. The 
corresponding events for distinct time-intervals are independent, whence the tail of T is 
no greater than geometric. □ 

The above recipe is exactly that of [21j when fi is monotone, and that of [L5\ when n 
is anti-monotone. 

Let G = {V,E) be a finite graph, and W ^ V a. non-empty set of vertices with \W\ 
even. Let r = (rg : e G -E) be a vector of numbers from (0,1], and let 0r,g be the 
random-cluster measure on G with edge-parameters r and q > 1. We write 0^ for 0r,g 
conditioned on the event that the open graph is W^-even, and note that (j)^^ is neither 
monotone nor anti- monotone. The event Sfj_ is easily seen to be increasing, and 1 G 5^. 
We may therefore apply Theorem 14.31 to the measure fi = (j)^^. 

Certain natural questions arise over the implementation of the above algorithm, and 
we shall not investigate these here. First, it is convenient to have a quick way to calculate 
a and (3 in (14.21) . A second problem is to determine the mean runtime of the algorithm, 
for which we remind the reader of the arguments of |2T| Sect. 5]. 

5 Random even subgraphs of planar lattices 

In this section, we consider random even subgraphs of the square and hexagonal lattices. 
We show that properties of the Ising models on these lattices imply properties of the 
random even graphs. In so doing, we shall review certain known properties of the Ising 
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model, and we include a 'modern' proof of the established fact that the Ising model on 
the square lattice has a unique Gibbs state at the critical point. 

Let G = (y,E) be a planar graph embedded in M?, with dual graph = (Vd,£^d)) 
and write Cd for the dual edge corresponding to the primal edge e E E. [See [12] for an 
account of planar duality in the context of the random-cluster model.] Let p e (0, |] and 
let uj E Q = {0, 1}^ have law 02p,2- There is a one-one correspondence between Q and 

= {0, 1}^'* given by uj{e) + ci;d(ed) = 1- It is well known that c^d has the law of the 
random-cluster model on G^ with parameters (1 — 2p) /{I —p) and 2, see [12] for example. 
For A <0 V, the boundary of A is given by dA = {v E A : v ^ w for some w ^ A}. [A 
similar notation was used in a different context in (11.21) . Both usages are standard, and 
no confusion will arise in this section.] 

For u E Q, let /o, /i, . . . , /c be the faces of (V^, r]{uj)), with /o the infinite face. These 
faces are in one-one correspondence with the clusters of (Vd, //(cUd)), which we thus denote 
by Kq, Ki, . . . , Kc, and the boundaries of the finite faces form a basis of £ = £{V,ri{uj)). 
More precisely, the boundary of each finite face /j consists of an 'outer boundary' and zero, 
one or several 'inner boundaries'; each of these parts is a cycle (and two parts may have 
up to one vertex in common). If we orient the outer boundary cycle counter-clockwise 
(positive) and the inner boundary cycles clockwise (negative), the face will always be on 
the left side along the boundaries, and the winding numbers of the boundary cycles sum 
up to 1 at every point inside the face and to outside the face. It is easy to see that the 
outer boundary cycles form a maximal family of independent cycles of {V,t]{uj)), and thus 
a basis of S; another basis is obtained by the complete boundaries Ci of the finite faces. 
We use the latter basis, and select a random subset of the basis by randomly assigning (by 
fair coin tosses) + and — to each cluster in the dual graph (Vd, 77(a;d)), or equivalently to 
each face fi of {V,r]{u)). We then select the boundaries Ci of the finite faces fi that have 
been given a sign different from the sign of the infinite face /q. The union (modulo 2) of 
the selected boundaries is by Proposition 12.21 and Theorem 13.11 a random even subgraph 
of G with parameter p. On the other hand, this union is exactly the dual boundary of the 
+ clusters of Gd, that is, the set of open edges e E E with the property that one endpoint 
of the corresponding dual edge Cd is labelled -|- and the other is labelled — . [Such an edge 
Cd is called a. +/— edge.] 

It is standard that the +/— configuration on Gd is distributed as the Ising model on 
Gd with parameter (3 satisfying 

^ = l-e-^ = ^^^^. (5.1) 

1 — p 1 + tanh p 

In summary, we have the following. 

Theorem 5.2. Let G be a finite planar graph with dual G^- A random even subgraph of 
G with parameter p E (0, |] is dual to the +/— edges of the Ising model on G^ with (3 
satisfying (15. ip . 

Much is known about the Ising model on finite subsets of two-dimensional lattices, and 
the above fact permits an analysis of random even subgraphs of their dual lattices. The 
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situation is much more interesting in the infinite-volume hmit, as follows. Let G = {V, E) 
be a finite subgraph of Z^, with boundary dV when viewed thus. A boundary condition 
on dV is a vector in Sgy = {— 1, +1}^^. For given 77 G S^y, one may consider the Ising 
measure, denoted Hy, on G conditioned to agree with r] on dV. We call any subsequential 
weak limit of the family {iiy : rj G S^y, \^ C Z^} a (weak limit) Gibbs state for the Ising 
model. It turns out that there exists a critical value of jS, denoted /3c, such that there is 
a unique Gibbs state when /3 < /3c, and more than one Gibbs state when /? > /?c. 

Consider the case when G is a box in the square lattice Z^. That is, G = Gm,n is the 
subgraph of Z^ induced by the vertex-set [—m, m] x [— n, n], where m,n & Z+ and [a, b] is 
to be interpreted as [a, 6] fl Z. It is a mild inconvenience that Gm,n is not an even graph, 
and we adjust the 'boundary' to rectify this. For definiteness, we consider the so-called 
'wired boundary condition' on Gm,n, which is to say that we consider the random- cluster 
measure on the graph „ obtained from Gm,n by identifying as one the set of vertices 
lying in its boundary dGm,n- 

It has been known since the work of Onsager that the Ising model on J? with parameter 
(3 is critical when e^^ = 1 + v^, or equivalently when the above random-cluster model on 
the dual lattice has parameter satisfying 

1-2P V2 



1 -p 1 + ^2 
that is, p = Pc where 

Pc = 7= = 1 - 5.3 

2 + V2 V2 

The Ising model has been studied extensively in the physics literature, and physicists 
have a detailed knowledge of the two-dimensional case particularly. There is a host of 
'exact calculations', rigorous proofs of which can present challenges to mathematicians, 
see [31 [20]. 

We shall use the established facts stated in the following theorem. The continuity of 
the magnetization at the critical point contributes to the proof that the re-scaled boundary 
of a large spin-cluster of the critical Ising model converges weakly to the Schramm-Lowner 
curve SLE3, see pl[25]. 

Theorem 5.4. The critical value of the Ising model on the square lattice is (3c = (3sd 
where /?sd = |log(l + v^) is the 'self-dual point'. The magnetization (and therefore the 
corresponding random- cluster percolation-probability also) is a continuous function of (3 
on [0, 00). 

We note the corollary that the wired and free random-cluster measures on J? are 
identical for p G [0, 1]; see [12^ Thms 5.33, 6.17]. 

Proof. These facts are 'classical' and have received much attention, see [20] for example; 
they may be proved as follows using 'modern' arguments. Recall first that the magnetiza- 
tion equals the percolation probability of the corresponding wired random-cluster model. 
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and the two-point correlation function of the Ising model equals the two-point connectiv- 
ity function of the random-cluster model (see [12] )• We have that P^d < Pc, by Theorem 
6.17(a) of [12] or otherwise, and similarly the random-cluster model with free boundary 
condition has percolation-probability whenever either (3 < f3sd oi (3 < (3c- 

By the results of [21 [22], the two-point correlation function Ti^laxCTy) of the spins at 
X and y decays exponentially as \x — y\ oo when j3 < (3c, and it follows by the final 
statement of [13] or otherwise that Pc = Ad- 

The continuity of the magnetization at /5 7^ /5c is standard, see for example [121 
Thms 5.16, 6.17(b)]. When /? = Pc, it suffices to show that the ± boundary-condition 
Gibbs states vr^^ and the free boundary-condition Gibbs state vr^^ satisfy vr^^ = vr^^ = vr^^. 
Suppose this does not hold, so that vr^^ 7^ vr^^ 7^ n^^. By the random-cluster representation 
or otherwise, the two-point correlation functions 7r'^^{ax(Jy) are bounded away from for 
all pairs x, y of vertices. By the main result of [U [16] (see also [10]) and the symmetry 
of TT^^, we have that n^^ = ^tt^^ + ^tt^^, whence 7ip^{axCry) is bounded away from 0. By 
[T2| Thm 5.17], this contradicts the above remark that the percolation-probability of the 
free-boundary condition random-cluster measure is at /5 = Psd = Pc- □ 

We consider now the so-called thermodynamic limit of the random even graph on ^ 
as m, — > 00. It is long established that the (free boundary condition) Ising measure 
on Gm,n converges weakly (in the product topology) to an infinite-volume limit measure 
denoted vr^. This may be seen as follows using the theory of the corresponding random- 
cluster model on (see [12] )• When P < Pc-, the existence of the limit follows more or 
less as discussed above, using the coupling with the random-cluster measure, and the fact 
that the percolation probability of the latter measure is whenever P < Pc- We write tt^ 
for the limit Ising measure as m, n ^ 00. 

The thermodynamic limit is slightly more subtle when P > Pc-, since the infinite- 
volume Ising model has a multiplicity of Gibbs states in this case. The (wired) random- 
cluster measure on „ converges to the wired limit measure. By the uniqueness of 
infinite-volume random-cluster measures, the limit Ising measure is obtained by allocating 
random spins to the clusters of the infinite-volume random-cluster model (see Section 4.6 
of [12])- Once again, we write Tijs for the ensuing measure on {—1, +1}^ , and we note that 

~ ~^ where vr^ denotes the infinite-volume Ising measure with ± boundary 
conditions. 

It has been shown in [B] (see also [HI Cor. 8.4]) that there exists (with strictly positive 
TT/j-probability) an infinite spin-cluster in the Ising model if and only if /? > Pc- More 
precisely: 

(a) if /? < Pc, there is vT/j-probability 1 that all spin-clusters are finite, 

(b) ii P > Pc, there is TT/j-probability 1 that there exists a unique infinite spin-cluster, 
which is equally likely to be a + cluster as a — cluster. Furthermore, by the main 
theorem of [8j or otherwise, for any given finite set 5* of vertices, the infinite spin- 
cluster contains, vr/j-a.s., a cycle containing S in its interior. 
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On passing to the dual graph, one finds that the random even subgraph of G^.^ with 
parameter p e (0, |] converges weakly as m, ri oo to a probability measure Pp that is 
concentrated on even subgraphs of and satisfies: 

(a') if p > Pc, there is pp-probability 1 that all faces of the graph are bounded, 

(b') if p < pc, there is pp-probability 1 that the graph is the vertex-disjoint union of 
finite clusters. 

(Note that (15.11) defines as a decreasing function of p, so the order relations are reversed.) 

We have thus obtained a description of the weak-limit measure Pp when p < |, and 
we note the phase transition at the parameter- value p = Pc- When p > |, a random even 
subgraph of „ is the complement of a random even subgraph with parameter 1—p. [It 
is a convenience at this point that „ is an even graph.] Hence the weak-limit measure 
Pp exists for all p G [0, 1] and gives meaning to the expression "a random even subgraph on 

with parameter p" . [It is easily verified that pi equals the measure defined in Theorem 
l2.6l for Z^, and thus describes a uniform random even subgraph of Z^.] There is a sense in 
which the random even subgraph on Z^ has two points of phase transition, corresponding 
to the values pc and 1 — Pc- 

We consider finally the question of the size of a typical face of the random even graph 
on Z^ when p^ < p < |. When p > Pc, this amounts to asking about the size of a 
(sub)critical Ising spin-cluster. Higuchi [T7] has proved an exponential upper bound for 
the radius of such a cluster, and this has been extended by van den Berg [4| to the 
cluster- volume. Thus, the law of the area of a typical face has an exponential tail. 

The picture is quite different when the square lattice is replaced by the hexagonal 
lattice H. Any even subgraph of EI has vertex degrees and/or 2, and thus comprises 
a vertex-disjoint union of cycles, doubly infinite paths, and isolated vertices. The (dual) 
Ising model inhabits the (Whitney) dual lattice of H, namely the triangular lattice T. 
Once again there exists a critical point pc = Pc(T) < | such that the random even 
subgraph of EI satisfies (a') and (b') above. In particular, the random even subgraph has 
a.s. only cycles and isolated vertices but no infinite paths. Recall that site percolation on 
T has critical value |. Therefore, for p = ^, the face of the random even subgraph 
containing the dual vertex x corresponds to a critical percolation cluster. It follows that 
its volume and radius have polynomially decaying tails, and that the boundary of F^, 
when conditioned to be increasingly large, approaches SLEg. See [221 [21] [SIE7|. The 
spin-clusters of the Ising model on T are 'critical' (in a certain sense described below) 
for all p G (pc(Tr), ^j, and this suggests the possibility that the boundary of F^, when 
conditioned to be increasingly large, approaches SLEg for any such p. This is supported 
by the belief in the physics community that the so-called universality class of the spin- 
clusters of the subcritical Ising model on T is the same as that of critical percolation. 

The 'criticality' of such Ising spin-clusters (mentioned above) may be obtained as 
follows. Note first that, since /? < Pc, there is a unique Gibbs state vr^ for the Ising model. 
Therefore, tt/j is invariant under the interchange of spin-values —1 +1. Let Rn be a 
rhombus of the lattice with side-lengths n and axes parallel to the horizontal and one of 
the diagonal lattice directions, and consider the event An that Rn is traversed from left 
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to right by a + path (i.e., a path v satisfying Oy = +1 for all y G z/). It is easily seen that 
the complement of An is the event that -R„ is crossed from top to bottom by a — path 
(see [TT], Lemma 11.21] for the analogous case of bond percolation on the square lattice). 
Therefore, 

7Tp{An) = i, 0<f3<(3c. (5.5) 
For X G Z^, let denote the spin-cluster containing x, and define 

rad(S'3;) = max{|2; — x\ : z E S}, 

where \y\ is the graph-theoretic distance from to y. By (15. 5p . there exists a vertex x 
such that 7r^(rad(S'a;) > n) > {2n)^^. By the translation-invariance of vr^, 

7r^(rad(S'o) >n) > i-, < /3 < (3,, (5.6) 

where denotes the origin of the lattice. The left side of (15. 6p tends to as n ^ oo, and 
the polynomial lower bound is an indication of the criticality of the model. 
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